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The critical behaviour of the Fully Frustrated XY model in presence 

of weak positional disorder is studied in a square lattice by Monte Carlo 

methods. The critical exponent associated to the divergence of the chiral 

correlation length is found to be v ~ 1.7 already at very small values of 

disorder. Furthermore the helicity modulus jump is found larger than the 

universal value expected in the XY model. 

The study of the fully frustrated XY (FFXY) modelBi in a planar lattices has attracted 
much attention since it provides a simple model for many physical systems. In particular it is 
related to Josephson-junction arrays and superconducting wire networks in a magnetic field! 
for which recent experiments show interesting properties!. In these systems a perpendicular 
magnetic field B penetrates the array inducing a quenched vortex density. When the field 
is such that each plaquette of area S in the planar lattice traps half a vortex, e.g. the ratio 
of the magnetic field to the flux unit / = BS/(p = 1/2, the system can be described by the 
FFXY model. As first noted by Villain! the FFXY model posses, besides the continous U(l) 
symmetry, a discrete Z 2 symmetry which is not present in the unfrutrated XY model. The 
existence of two different symmetries leads to the possibility of two critical temperatures that 
have been investigated by both analytical methodsHi and Monte Carlo (MC) simulations!~0. 
However, in spite of accurate MC simulations, the scenario is not fully understood. In 
fact, there are numerical evidences supporting both the existence of two very close critical 
temperatures (T^ 2 > T^ 1 -*) with critical behavior typical of Ising and Thouless-Kosterlitz- 
Berezinski (TKB) transitions respectively^^! and the existence of a single transition with 
novel critical behaviorB. In the latter case it has been claimed that the transition associated 



1 



to the Z2 symmetry is characterized by a correlation length critical exponent v lower than 
that expected in the Ising case. 

Recent experiments on Nb wire networks^ in an external magnetic field have shown 
that the critical exponent v which can be extracted by the scaling of current-voltage (I-V) 
characteristics, is in agreament neither with a Ising-like exponent iy = 1) nor with the value 
suggested in Ref.P ] [y ~ 0.83). On the contrary, the scaling of the experimental data 
suggests a much larger value. 

The aim of the present work is to investigate if the effect of weak quenched disorder in 
the FFXY model is able to change the critical behaviour of the system and to explain the 
unexpected critical behaviour. 

The hamiltonian of the FFXY model is 

H = -Jj2cos(6 i -9 j -A id ), (1) 

(ij) 

where Q{ is the phase in the site i and 
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A,j = — A-dl (2) 

is the line integral of the vector potential along the line (z — j) and 4>q is the flux quantum. 
In order to simulate positional disorder we consider a uniform distribution of the Aij whose 
average value (Aij) is equal to 7r every second column and zero otherwise. The disorder 
configuration is thus specified by the probability distribution 

p(Ati) \± -*.<(*,-<**»<*. (3) 

otherwise 

The parameter m controls the degree of randomness and is related to the variance according 
to a 2 = 

Of course, the presence of positional disorder effects the behaviour of the model. In fact 
it has been shown by renormalization group techniques^ and Monte Carlo simulationsEi that 
the critical temperatures associated to both Z 2 and U(l) simmetries decrease with increasing 
disorder. It has been also suggested for the U(l) transition a reentrant behaviour, even if 



recent work on the disordered unfrustrated XY model seems to exclude this possibility^. 
However, the critical behaviour that characterize the Ising-like and the TKB transitions 
have not been studied to our knowledge. In particular, the study of the critical exponent v 
associated to the divergence of the coherence lenght of the Ising-like transitions is interesting 
in view of its possible connection with the I — V characteristics reported by LingS. In fact 
the experimental data I-V, or equivalently the data J-E, scale very well according to the 
relation^ 

E = J | T-T c \ uz G(J | T-T c \~ v ) (4) 

where the dynamical critical exponent z is found to be about 2 and v = 1.7 ± 0.5. This 
type of scaling is consistent with the reasonable assumption that the relevant excitations 
that contribute to the critical behaviour of the I-V characteristic are domain-wall. 

In order to obtain an estimation of v we have calculated the Binder's cumulant of the 
staggered magnetization M 

TJ 1 < M4 > ™ 
where M, associated to the chirality of each plaquette, is defined as 

M = j2 1 EMr^- 1 • ( 6 ) 

r 

In eq. @ the chirality m r is the rotation of the current around a plaquette 

m = —=(sin<f)i2 + sm023 + sm0 34 + sm0 41 ) (7) 
V8 

where 4>ij is the phase difference at the edges of the plaquette. 

Binder's cumulant allows an accurate estimation of the critical temperature associated 
to the ordering of the plaquette chirality (Ising-like transition). In fact U(T C ,L) does not 
depend on lattice size L for large systems and, then, T c can be identified without making any 
assumption on the critical exponents. Once a satisfactory estimation of T c is obtained the 
critical exponent v is estimated through a data collapsing with the parameter v only. The 



Binder's cumulant U has been calculated by Monte Carlo (MC) simulations using a standard 
Metropolis dynamics. For each disorder configuration the mean value (M) has been obtained 
averaging over 4 x 10 5 MC sweeps after discarding the first 10 5 for equilibrating the system. 
Then the mean values are further averaged over a number of disorder configurations that 
range from 200 to 400 depending on lattice size and amount of disorder present in the system. 

In Fig. 1) it is shown how the critical temperature T^ 2 can be extracted by using the MC 
estimations of U for different lattice sizes. In the case reported in the figure, if we exclude 
the smallest lattice size, L = 12, all the data cross about Tf 2 = 0.438 ± 0.003. From this 
type of analysis it is, then, possible to follow T^ 2 as a function of the disorder (Fig. 2)). The 
results confirm that T^ 2 decrease with increasing disorder going from T^ 2 = 0.452 ± 0.002 
for the case without disorder (0 m = 0) to T^ 2 = 0.350 ± 0.005 for the case with largest 
disorder studied (0 m = 0.245). The results are in qualitative agreement with those reported 
by Choi et al.0. For larger values of the disorder preliminar results seems to suggest a 
continous collaps of T^ 2 . In Figs. 3a) and 3b) we show Binder's cumulant U as a function 
of | T — T C Z2 | L l l u for two values of m . By using the estimated values of T^ 2 at m = 0.098 
and m = 0.196 we see that the best data collapsing is obtained for v = 0.9 ± 0.2 in the 
case of lower disorder and for v = 1.7±0.2 in the other case. Even if the numerical errors 
do not allow a very accurate estimation of v it is very clear that already for very weak 
disorder (<f) m > 0.19) the critical exponent is quite different from that expected in the Ising 
model [y = 1) and from that suggested by Granato et al.i in the case without disorder. 
Instead the value obtained is very close to that extracted from the I — V characteristics!. 
It is also interesting to note that a very similar value is obtained in the 2d XY spin glass 
that has a critical temperature T^ 2 = 0.0^0 The same analysis for m = and m = 0.245 
gives clearly v ~ 1 and v ~ 1.7, respectively. On the other hand for the intermediate value 
<f) m = 0.147 the data are not able to clearly exclude either v ~ 1 or v ~ 1.7; larger lattices 
and more MC sweeps are needed to reach a satisfactory estimation of v. These results 
indicate the existence of a threshold disorder above which v changes from 1 to 1.7. The 
value of the threshold can be, then, estimated to be m ~ 0.19 but it must be understood 



as an upper limit. 

As in the FFXY model we expect the system to present a jump in the helicity modulus 
at temperature lower or equal to T/ 2 . In order to study this possibility we have calculated 
the helicity modulus that for the system studied takes the form 

r = ^ j-^<(l>(^ - Q 3 ~ A ^ ■ ) 

+ (H cos(0i - 9j - Aij) {tij ■ xf) 

m 

+ ^(YtMOi-Oj-A^ij-zA (8) 

where e^- and x are the unitary vectors in the directions of the link % i— > j and x, respectively. 
The symbol (•) indicates the average over the spin configurations. The elicity modulus V 
has to be averaged further over the disorder configurations to obtain the mean value F. 
Following the analysis developed for the XY mo deE and extended to the FFXY modelEI 
we can estimate the critical temperature T^ 1 ' by using the following ansatz for the size 
dependence of T 

^=T + - (9) 

2T C 2(lnL + ln/ ) K) 

where To = 1 corresponds to the universal jump expected at the critical temperature in 

a KTB transition. This critical scaling, being based on the mapping between a neutral 

Coulomb gas and a XY model, is valid only if the transition is of KTB type. Therefore eq. 

(§) can be used as a test for the existence of a standard KTB transition. 

In Figs. 4a) and 4b) we show the results of this analysis for two different values of the 

disorder. In both cases it is not possible to obtain a satisfactory fit of the data by using 

eq. (|9|) with T = 1. However, a more accurate data analysis shows that a satisfying fitting 

of the data can be obtained assuming a helicity modulus jump larger than the universal 

value. From such an analysis we get TfW = 0.43 ± 0.003 and T = 1.2 for <\> m = 0.098 

and T^ 1 ) = 0.38 ± 0.003 and T = 1.7 for <p m = 0.196. In both cases the fitting parameter 

Iq takes the same value found in the FFXY model without disorder, l = 0.27,fHl! and T^ 1 ) is 



slightly smaller than T^ 2 even if the estimated numerical error does not rule out T^W = Tf 2 
in the case of lowest disorder (0 m = 0.098). ill This analysis shows that the only way to fit 
the data with eq.@ is to assume a larger helicity modulus jump and this is an indication 
that the transition associated to this jump is not a KTB transition. 

Summarizing we have calculated, by MC simulations, the critical temperature and the 
critical exponent v associated to the chiral degree of freedom in a disordered FFXY model. 
Unexpected we find that the critical exponent v already for very weak disorder becomes 
equal to the value found in the XY spin glass, i.e. v ~ 1.7. Correspondently, the helicity 
modulus jump becomes larger than the universal value expected in KTB transition. Finally, 
we note that the value obtained for v is in agreement with recent results in Josephson wire 
networks .0 
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FIGURE CAPTIONS 



Fig. 1 The Binder' s cumulant U (see eq.([5])) vs. temperature T for different sizes L. At 
the critical temperature U does not depend on L for large L. 

Fig. 2 The critical temperature T^ 2 as a function of the parameter <p rn that controls the 
degree of disorder present in the system. The continous cirve is a guide for the eye 
only. 

Fig. 3 The Binder's cumulant (see eq.(|J)) as a function of (T — T ( f 2 )L ( ^ 1//iy - ) for two values of 
the parameter that controls the disorder <p m . In Fig. 3a) the data for m = 0.098 are 
reported having choosen T^ 2 = 0.438 and v — 1. In Fig.3b) the data for m = 0.196 
are reported having choosen T^ 2 = 0.395 and v = 1.7. Temperature is in units of 
k B /J. 

Fig. 4 The scaled helicity modulus ^f- as a function of InL for two values of the pa- 
rameter that controls the disorder <p m . In Fig. 4a) the data for cf) m = 0.098 and 
T = 0.45, 0.44, 0.43, 0.42 (from below to the top) are reported. The best fit is given 
by the full curve that corresponds to eq.([5]) with r = 1.2 and Iq = 0.27. In Fig. 4b) 
the data for m = 0.196 and T = 41, 39, 38, 37 (from below to the top) are reported. 
The best fit is given by the full curve that corresponds to eq.(|9|) with r = 1.7 and 
L = 0.27. 
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